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Abstract 

We describe the use of doubly-unintegrated parton distributions in hadron-hadron 
collisions, using the (z, /^-factorisation prescription where the transverse momentum of 
the incoming parton is generated in the last evolution step. We apply this formalism to 
calculate the transverse momentum (Pt) distributions of produced W and Z bosons and 
compare the predictions to Tevatron Run 1 data. We find that the observed Pt distri- 
butions can be generated almost entirely by the leading order q\ q 2 — > W, Z subprocesses, 
using known and universal doubly-unintegrated quark distributions. We also calculate 
the Pt distribution of the Standard Model Higgs boson at the LHC, where the dominant 
production mechanism is by gluon-gluon fusion. 



1 Introduction 



At present, it is not straightforward to describe the transverse momentum (Pt) distributions of 
electroweak bosons produced in hadron-hadron collisions. In the usual collinear approximation, 
the transverse momentum of the incoming partons is neglected and so, for the Born level 
subprocesses q\ q 2 — > V (where V = 7*, W, Z) or g 1 g 2 — > H, the transverse momentum of the 
final electroweak boson is zero. Therefore, initial-state QCD radiation is necessary to generate 
the Pt distributions. Both the leading order (LO) and next-to-leading order (NLO) differential 
cross sections diverge for Pt <C Mv,h, with terms proportional to \n(Mv,H / 'Pt) appearing due 
to soft and collinear gluon emission, requiring resummation to achieve a finite Pt distribution. 

Traditional calculations combine fixed-order perturbation theory at high Pt with either 
analytic resummation or numerical parton shower formalisms at low Pt, with some matching 
criterion to decide when to switch between the two. In addition, a parameterisation is needed 
to account for non-perturbative effects at the lowest Pt values. Analytic resummation can be 
performed either in the transverse momentum space (see, for example, Ij) or in the Fourier 
conjugate impact parameter space (see, for example, 0). 



An alternate description is provided in terms of unintegrated parton distribution functions 
(UPDFs), where each incoming parton carries its own transverse momentum k t , so that the 
subprocesses qi q 2 — > V and g x g 2 — > H already generate the LO Pt distributions in the so-called 
k t - factorisation approach. 1 It has been shown in [3] that UPDFs obtained from an approximate 
solution of the CCFM evolution equation jS] embody the conventional soft gluon resummation 
formulae. 

The UPDFs that we use are obtained from the familiar DGLAP-evolved PDFs determined 
from a global parton analysis of deep-inelastic and related hard-scattering data. The transverse 
momentum of the parton is generated entirely in the last evolution step jHllZj- Angular-ordering 
constraints are imposed which regulate the singularities arising from soft gluon emission, while 
the virtual terms in the DGLAP equation are resummed into Sudakov form factors. 2 In jH] 
it was shown that it is necessary to extend the 'last-step' formalism of 7\ to consider doubly- 
unintegrated parton distribution functions (DUPDFs) in order to preserve the exact kinematics. 
It was demonstrated that the main features of conventional higher order calculations can be 
accounted for within a much simpler theoretical framework, named (z, /cj-factorisation. 3 

Strictly speaking, the integrated PDFs used as input to the last evolution step should them- 
selves be determined from a new global fit to data using the (z, /^-factorisation approach. For 
the present paper, we take the input PDFs from a global fit to data using the conventional 
collinear approximation ^JJ. This treatment is adequate for these initial investigations. How- 
ever, we expect it to lower our predictions for quark-initiated processes by ~10% compared to 
the case where the input PDFs are determined from a global fit using the (z, fc^-factorisation 
approach. We will illustrate this point in Section 5 by comparing predictions for the proton 
structure function F 2 in the collinear approximation and in the (z, /^-factorisation formalism. 

The 'last-step' prescription has some features in common with the initial-state parton shower 
algorithms implemented in Monte Carlo event generators (for a recent review, see [T2j) such as 
the DGLAP-based herwig [J3] and pythia [JJ] programs and the CCFM-based cascade [T3] 
program. The main advantage of our approach is that we use simple analytic formulae which 
implement the crucial physics in a transparent way, without the additional details or tuning 
which are frequently introduced in Monte Carlo programs. The Pj> distributions are generated 
entirely from known and universal DUPDFs. For example, fits to Z production data at the 
Tevatron using parton showers favour a large intrinsic partonic transverse momentum (kt) ~ 2 
GeV, while confinement of partons inside the proton would imply a (kt) ~ 0.3 GeV [To] . 

The DGLAP-based parton showers used in El are theoretically justified only in the 
limit of strongly-ordered transverse momentum, since only the collinear divergent part of the 
squared matrix element is kept in each parton branching. Similarly, the CCFM-based parton 

^or an introduction to fc t -factorisation, see 

2 It was found in [7] that the 'last-step' prescription gave similar results whether the input PDFs were evolved 
with the DGLAP equation or with a unified BFKL-DGLAP equation |5J|, indicating that angular ordering is 
more important than small- a; effects. 

3 Hcre, z is the splitting fraction associated with the last evolution step, where a parton with (light-cone) 
momentum fraction x/z splits to a final parton with (light-cone) momentum fraction x. The idea is an extension 
of the original DDT formula 10 ; however, in comparison with |1U| we go beyond the double leading logarithmic 
approximation (DLLA) and account for the precise kinematics of the two incoming partons, as well as the 
angular ordering of emitted gluons. 



2 



shower used in jT3] is strictly justified only in the limit of strongly-ordered angles, which reduces 
to the limit of strongly-ordered transverse momentum as long as x is not too small. In this 
limit, the transverse momentum generated in all evolution steps prior to the last is negligible. 
Therefore, neglecting transverse momentum in every evolution step prior to the last should be 
a good approximation to the parton shower algorithms in which finite transverse momentum is 
generated at every evolution step. 

The goal of the present paper is to demonstrate that the Pt distributions of electroweak 
bosons can be successfully generated by DUPDFs. We do not aim to produce a better de- 
scription of the data than existing calculations, but rather a simpler analytic description which 
reproduces the main features. With this approach, it is easy to see the physical origin of the 
Pt distributions and to identify the most important Feynman diagrams. Since the DUPDFs 
are universal — that is, they apply equally well to all hard hadronic processes — it is important 
to check them in a new kinematic domain. 

One topical application is the prediction of the cross section for diffractive Higgs boson 
production at the LHC [T7], which is driven by the unintegrated 4 gluon distribution f g (x, hf, /x 2 ), 
where // is the hard scale of the subprocess. At the moment, the only possibility to check the 
behaviour of UPDFs in the domain k t <C /i is to compare predictions with the observed Pt 
distributions of W and Z bosons produced at the Tevatron. We will show that the doubly- 
unintegrated quark distributions, generated directly from the known integrated PDFs under 
the 'last-step' prescription, satisfactorily describe these data, including the region of interest, 
P T < M w ,z- 

In Section El we describe the formalism for (z, /^-factorisation at hadron-hadron colliders, 
and in Section El we apply it to calculate the Pt distributions of electroweak bosons. Numerical 
results are given in Sectional 

2 (z, /^-factorisation at hadron-hadron colliders 

We now extend the formalism of which concerned deep- inelastic scattering (DIS), to hadron- 
hadron collisions. The basic idea is illustrated in Fig.^a), which shows only one of the possible 
configurations. All permutations of quarks and gluons must be included. The arrows show 
the direction of the labelled momenta. The blobs represent the familiar integrated PDFs. 
The transverse momenta of the two incoming partons to the subprocess, represented by the 
rectangles labelled a qi q2 in Fig. Q are generated by a single parton emission in the last evolution 
step. 

We use a Sudakov decomposition of the momenta of the two incoming partons: 

k-i — Xi Pi Pi Pj + (1) 

where = (1, 2) or (2, 1). We work in the centre-of-mass frame of the colliding hadrons and 
neglect the hadron masses so that the squared centre-of-mass energy is s = (P1+P2) 2 — 2 P\-Pi- 

4 To be precise, the skewed unintegrated gluon distribution is required. However, in the relevant small- a; 
domain the skewed effect can be included by the Shuvaev prescription |18j . 
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(a) (b) 




Figure 1: (a) The transverse momentum of each parton entering the subprocess is generated 
by a single parton emission in the last evolution step, (b) Illustration of (z, A^)-factorisation: 
the last evolution step is factorised into f gi (xi, z iy kf t , where i — 1, 2. 
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Then, 5 

Pi = (P+, Pf, Pi ± ) = (1, 0, 0), P 2 = (0, 1, 0) A; i± = (0, 0, fe i)t ). (2) 

The penultimate propagators in the evolution ladder have momenta U = XiPi/zi, so that the 
partons emitted in the last step have momenta 

Pi = k - hi = — (1 - Zi) Pi + f3i Pj - k i± . (3) 
The on-shell condition for the emitted partons, p 2 = 0, determines 

ft - mt^) • (4) 

where = k 2 t /s, so that the two incoming partons have virtuality kf = —k 2 t /{l — Zi). The 
total momentum going into the subprocess labelled a qi 92 in Fig. ^ is 

q = ki + k 2 = (x 1 - p 2 ) Pi + (x 2 - Pi) P 2 + q±, (5) 

where q± = ki± + k 2 ±. The kinematic variables obey the ordering 

< Pj < Xi < Zi < 1. (6) 

If pi are gluon momenta, then we must additionally impose angular-ordering 6 constraints 
due to colour coherence [5]: 

6 < s < 6, (7) 

where & = p^ /pf and H = q~ /q + ■ That is, the subprocess separates gluons emitted from each 
of the two hadrons. This condition leads to a suppression of soft gluon emission: 

Zi < (8) 



with fix = XiVsE and /i 2 = x 2 y s / E. 

The partonic cross sections, represented by the rectangles in Fig. HI depend on Xi, Zi, and 
kij through the momenta of the incoming partons k,i ([TJ. Therefore, the partonic cross sections 
must be convoluted with parton distributions depending on the same variables. We can define 
DUPDFs f a (x, z,k 2 , n 2 ) which satisfy the normalisation conditions 



dz I -jjffa{x,z, k 2 ,n 2 ) = a(x,/i 2 ), (9) 



t 

2\ _ „ „(„ ,.2\ m „f„ ,,2\ 



for fixed x and /i, independent of z and k t . Here, a(x,/i) = xg(x,pr) or xq(x,fx) are the 
conventional integrated PDFs. Explicit formulae for the DUPDFs are given in The doubly- 
unintegrated quark distribution is 



t I J 2 2\ rn 11 2 2\ a s{kt ) 

f q {x, z,k t ,fi)= T q (k t , fi ) — 



(10) 



5 The 'plus' and 'minus' light-cone components of a 4-momentum P are P ± = P° ± P 3 . 
6 "Angular ordering" is in fact a misnomer. It is rapidity ordering which should be applied. 



where P qq (z) and P qg (z) are the unregulated LO DGLAP splitting kernels, and the quark 
Sudakov form factor is 



T g ( k hv 2 ) = ex P 



'&k 2 a s (K 2 t ] 



2tt 



d£P«(C) • 



(11) 



The doubly-unintegrated gluon distribution is 



f g (x, z, k 2 , fi 2 ) — T g (k 2 , fi 2 



2vr 



J2P 9q (z)- zq (- z ,k 2 ) + P 99 (z)* z g (?-,k 2 )Q 



fji + kt 



where the gluon Sudakov form factor is 



T g( k t^ 2 ) = ex P 



3.K 2 a s (K 2 ] 



2tt 



dCC^(C) + ^ / d(P qg (() 



— z 



(12) 



(13) 



with rip the number of active flavours. 

In terms of the DUPDFs, the hadronic cross section a is related to the partonic cross sections 
a ai a2 by the (z, /c 4 )-factorisation formula 



0,1,0,2 



1 dx\ f l dx 



o X\ J x 2 



dzi / dz< 



X2 



°° dk 2 



>dkl 

p I 7.2 

"a,* JO ^2,7; 



/ai (^lj ^l, ki +, Hi) fa 2 (x 2 ,z 2 ,k 2t ,/i 2 ) a 



a± 02 



(14) 

This formula is represented schematically in Fig. E^b) for the case where the partons a\ and a 2 
are both quarks. The partonic cross sections in (JT4*j) are given by 



da 



Ol C12 



d$ aia2 |.M aia2 | / F ai 



. 2 



(15) 



where the flux factor F aia2 = 2x\X 2 s. The last evolution steps in Fig. ^a) factorise from 
the rest of the diagram, to give the LO DGLAP splitting kernels, in the leading logarithmic 
approximation (LLA) where only the dkf t /k 2 t term is kept. Therefore, |.M aia2 | 2 is calculated 
with the replacement ki — > Xi Pi in the numerator in order to keep only this collinear divergent 
term. However, any propagator virtualities appearing in the denominator of |.M aia2 | 2 may be 
evaluated with the full kinematics, as may the phase space element d$ ai a2 . 

For this approach to work, it is vital that the dkf t /kf t term is obtained only from ladder- 
type diagrams like that in Fig. EJa), and not from interference (non-ladder) diagrams. This 
is true if we use a physical gauge for the gluon, where only the two transverse polarisations 
propagate. For hadron-hadron collisions, the natural choice is the planar gauge where the sum 
over gluon polarisations is performed using 



"txiA'v = ~9uv H j : 

k ■ n 



(16) 



where we take the gauge-fixing vector n = x\ P\ + x 2 P 2 . Such a gauge choice ensures that the 
dk 2 t /k 2 t term is obtained from ladder-type diagrams on both sides of the subprocess represented 
by the rectangle in Fig. Ufa) [TU] . 
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A related requirement is that terms beyond the leading dk 2 t /k 2 t term, coming from non- 
ladder diagrams, for example, give a negligible contribution. Such terms are proportional to 
the Sudakov variable /3j (@J) and hence vanish in the limit that Zj — > or kf t /s — > (unless 
Zi — > 1). Away from these limits it is not obvious that these "beyond LLA" terms will be small, 
a necessary condition for the factorisation to hold. For the case of inclusive jet production in 
DIS and working in an axial gluon gauge, it was observed in [§| that the main effect of the 
extra terms was to suppress soft gluon emission. When the angular-ordering constraint (jHJ) was 
applied, the extra terms were found to make a negligible difference to the cross section. For 
hadron-hadron collisions, although the number of possible non-ladder diagrams is larger, it is 
therefore reasonable to expect that the extra terms will have little numerical effect, at least 
for ki t t less than the hard scale of the subprocess. A similar argument is made to justify the 
approximation made in the DGLAP-based parton showers used in Monte Carlo simulations, 
where only the collinear divergent part of the squared matrix element for each parton branching 
is kept and angular ordering is imposed in all evolution steps to account for some of the missing 
terms. Here, we are more conservative and apply this approximation to the last evolution step 
only. 

The DUPDFs in (|14|) are only defined for k itt > fi , where Ho ~ 1 GeV is the minimum 
scale for which DGLAP evolution of the integrated PDFs is valid. The approximation of the 
ki,t < Ho contribution made in [§| was to take the limit k^t — > in the kinematic variables (and 
in a aia2 ), then to make the replacement 

rl f^odk 2 ]- 

/ C ^ Zi / 7,2 fa t { X ii Z ii^i,ti Hi) = a i { X i i Ho) T ai {Hoi Hi ) , (17) 
Jxi JO K i,t 

where T ai (fil, fi 2 ) are the Sudakov form factors ((TTj) or (fTHjl . This replacement ensures that 
the normalisation conditions © are satisfied (see [9J. A better approximation, which retains 
the ki >t dependence, is to take the limit Z{ — > in the kinematic variables, then to make the 
replacement 

/•: _ //^ 

/ dZj f ai (Xi, Zi, k i ^ t , Hi) = faX X ii^ipHi) = 2 ai \ Xi i ^o) -^Oj (Hoi Hi ) ) (1§) 
Jxi Ho 

where we have used (23) of 0. The requirement that f ai (xi,kj t ,fi 2 ) ~ kf t as k\ t — ► is a 
consequence of gauge invariance (THj. 

A more complicated extrapolation of the DUPDFs for k t < fi , which allows both the k t 
and z dependence to be retained in the kinematic variables, is to assume the polynomial form 



2 2 kt 
fa(x, Z, k t , H ) = ~~2 

Ho 



k 2 

A(x, z, h 2 ) + -\ B(x, z, h 2 ) 
l^o 



(19) 



The two coefficients A and B can be determined to ensure continuity at k t = fi and the correct 
normalisation (JHj), leading to 

A(x, z, /i 2 ) = -f a (x, z, hI, H 2 ) + 2 a& Ho) T a(Hoi ^ 2 ) I i 1 ~ x )i ( 20 ) 
B(x, z, /j, 2 ) = 2 f a (x, z, hI, H 2 ) - 2 a(x, hI) T a(Hai H 2 ) I (! - x )- ( 21 ) 
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Figure 2: LO Feynman diagram contributing to the Pt distribution of W or Z bosons in the 
(z, /^-factorisation approach. 




H 



q = ki + k 2 



Figure 3: LO Feynman diagram contributing to the Pt distribution of Higgs bosons in the 
(z, k t )-factorisation approach. The cross represents the effective vertex in the limit Mu <C 2m t . 

If fjTHj) becomes negative, which is possible for k t <C fi®, then f a (x, z,k 2 , /z 2 ) is simply set to 
zero. 

This extrapolation of the perturbative formulae accounts for some non-perturbative "in- 
trinsic" kt of the initial partons, which is often parameterised by a Gaussian distribution with 
(k t ) < (j,q. Numerical results are insensitive to the precise form (JT7J), (jTSjl . or (jT^jl used for the 
k t < /i contribution. However, we will use the form (JT§j) which ensures continuity at k t = fi . 



3 Application to the Pt distribution of electroweak bosons 

Perhaps the simplest application of (z, fc 4 )-factorisation at hadron-hadron colliders is elec- 
troweak boson production, where at LO the subprocess is simply g* q% — > V, illustrated in 
Fig. |21 or g\g\ — > H, illustrated in Fig. 01 Here, the '*' indicates that the incoming partons, 
with momenta fc, given by ((TJ, are off-shell with virtuality kf = —kf t /(l — Zj). In the collinear 
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approximation, these diagrams give the Born level estimate of the total cross section a. When 
each parton carries finite transverse momentum fcj jt , the final electroweak boson has transverse 
momentum q t = k l t + fc 2 ,t, so we can calculate the Pt distribution 

{[n aS(q t -P T ), (22) 



with a given by the (z, /^-factorisation formula (fTf|) . 

The partonic differential cross sections for W or Z production are 

<!(T q*2 ^V) = ^V2G F M 2 V V 2 5{q 2 - M v ) 5{q t - P T ), (23) 



dP T ™ ^ ' N c 

where Nc = 3 is the number of colours, Gf is the Fermi coupling constant, Vyy = |lg ig2 | 2 is 
the CKM matrix element squared, and Vf = V 2 + A 2 is the sum of the vector and axial vector 
couplings squared. 

The dominant mechanism for SM Higgs production in hadron-hadron collisions is by gluon- 
gluon fusion via a top quark loop. For the case where Mh <C 2m*, the well-known effective 
ggH vertex can be derived from the Lagrangian (2H| 

where v 2 = (y/2 Gf)^ 1 , is the gluon field strength tensor, and H is the Higgs field. The 
partonic differential cross section for SM Higgs production is then 

^{9l9l - H) = ^ G Z* lX2S a 2 (M 2 H ) 5(q 2 - M 2 H ) 5(q t - P T ). (25) 
(xPt 57o 7r 

For q t = 0, these expressions (|23|) and (|25|) are exactly as in the collinear approximation. 
The difference arises when we consider the precise kinematics 

q 2 = s [( Xl - /3 2 )(x 2 -/3J-R], R = q 2 /s, (26) 
Qt = \ki,t + fc 2 ,t| 2 = k 2 1<t + kl t + 2 k ht k 2 , t cos0. (27) 

Applying the (z, fcf)-factorisation formula (fTlj) . the first delta function in (|2l?|) and (f2~K|) can 
be used to do the x 2 integration in (|14j) . while the second delta function can be used to do the 
k 2yt integration. In addition, we need to average over the azimuthal angle between fci jt and 
k 2 ,t- 

The final hadronic differential cross sections for W or Z production are 
<LPr N c ^ + ,^ + Jo *i J X1 J X2 J k 2 J 2tc k 2jt \k 2it + k 1:t cos^>| 



X2=X^ k2 t=fc ± 



2.1 



(xi, Zi, kl t , nl) f q2 (x 2 ,z 2 ,klt,fil), (28) 



x — 

xi x 2 - P\ f3 2 
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with t = My I 's, k 2t = —ki t t cos</> ±WP^- ^i,tsin 0, and 



_L T + fl + IlA + I? L 



± 4/ T + + + 



l-*2 



Z 2 T 2 

l-z 2 



(29) 



In practice, the kinematic constraints © mean that the x 2 = x 2 solution does not contribute. 
The corresponding result for SM Higgs production is 



da 



dP 



T 



^4(mde E 



«o &i J xx J*? 2 Jo k !,t Jo 2tt k 2 ,t\k 2 ,t + h,t cos, 

X2=x 2 k2,t=k 2 t 



X 



x x x 2 



f g (x 1 , zx, kl t , nl) f g {x 2 , z 2 , k 2t , (30) 



xx x 2 - fix (3 2 

where r = Mjj/s and k 2t and above. Note that we have taken the ggH vertex in 

the Mh <C 2m t limit. For Mh < 2m t , the correction to the total cross section due to the top 
quark mass can be approximated [21] by a factor 

on 2 



1 + 



M. 



H 



2m t 



(31) 



The ki t t < fio contributions of (J25j) and (jSOjl are accounted for using the approximation (fT9"j) . 7 



4 The if-factor 

In the collinear approximation, higher order QCD corrections to the LO diagrams, qx q 2 — > V 
or Pi P2 ~~ ► H, are known to be significant when calculating the total cross section. The ratio 
of the corrected result to the leading order result is the so-called .fT-factor. A part of these 
higher order corrections is kinematic in nature, arising from real parton emission, which we have 
already accounted for at LO in the (z, k t ) -factorisation approach (see Fig. [TJ. Another part 
comes from the logarithmic loop corrections which have already been included in the Sudakov 
form factors (jllj) and (|13j) . However, we need to include the non-logarithmic loop corrections 
arising, for example, from the gluon vertex correction to Figs. 121 and El 

A large part of these non-logarithmic corrections have a semi-classical nature and may be 
obtained from the analytic continuation of the double logarithm in the Sudakov form factors in 
going from spacelike (DIS) to timelike (Drell-Yan) electroweak boson momenta j22]- Accounting 
for the running of the strong coupling as, the Sudakov form factors in the DLL A are 

T^/i 2 ) =exp f-^L In (32) 

and 

T s (fc 4 V)=exp(-^LlnL), (33) 

7 Other contributions in the region kij < fJ.o, such as the inclusion of additional intrinsic partonic transverse 
momentum with (kt) ~ 0.3 GeV, would only affect the Pt distributions at very low Pt < 1 GeV. 
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where L = ln(/x 2 /AQ CD ), b = (33 — 2n F )/ (12ti), and the colour factors are Cp = 4/3 and 
C A = 3. 

Replacing /x 2 by — /x 2 , we obtain the 7r 2 -enhanced part of the i^-factors: 



and 



K(q* iq *^V) 



K(g* ig *^H) 





( 

~ exp 1 


T q (k 2 ^ 2 ) 


T 9 (kl-^ 2 ) 


( 

~ exp 1 


T 9 (k 2 ^ 2 ) 



2tt 



«s(/j 2 
2tt 



7T 



7T 



(34) 



(35) 



A particular scale choice /x 2 = P^ 3 My^ has been found (2H! to eliminate certain sub-leading 
logarithms in the Sudakov form factors. Therefore, we choose this scale to evaluate «s(/x 2 ) in 
flMD and dSHD- 



5 Numerical results 

5.1 W and Z boson production at the Tevatron 

The Pt distributions of produced W and Z bosons were measured by the CDF [21] and D0 
[2SH2H1 Collaborations during the Tevatron Run 1, in pp collisions at a centre-of-mass energy 
of y/s = 1.8 TeV. Measurements were made of W — > ev and Z — > ee decays; therefore, we 
must multiply the theoretical predictions for W or Z production by the appropriate leptonic 
branching ratios. 8 We use the MRST2001 LO PDFs [11J as input, with /xjj = 1.25 GeV 2 . The 
LO predictions for the Pt distributions (|28|). integrated over bins of 1 GeV, are shown by the 
dashed lines in Figs. [HEl and El The integrated luminosity uncertainty (3.9% for CDF or 4.4% 
for D0) is not included in the error bars for the plotted data. The LO predictions multiplied 
by the .fT-factor (jH^j) are shown by the solid lines. Although the .fT-factor makes up the major 
part of the discrepancy, we see that the solid lines in Figs. [H El and El still underestimate the 
precise measurements at small Pt to some extent. Normalising to the total measured cross 
sections (248 pb, 221 pb, and 2310 pb) by factors 1.23, 1.10, and 1.08 respectively, as shown by 
the dotted lines in Figs. [H El and El gives a very good description of the data over the entire 
Pt range. 

Although it makes sense to take /i 2 , as low as possible, the insensitivity of the Pt distributions 
to the k t < fio treatment can be demonstrated by taking /Xq = 2.5 GeV 2 . The Pt distributions 
obtained are practically identical to those with /Xq = 1.25 GeV 2 , with less than a 0.5% change 
in the total cross sections. 

Notice that the predicted Pt distribution of Z bosons peaks about 0.5-1.0 GeV below the 
CDF data (Fig. 0]). One possible explanation for this is provided by non-perturbative power 
corrections, part of which may be interpreted as a negative correction of about —3 GeV 2 to 
the factorisation scale at which the integrated PDFs are evaluated [2H]- Such a shift in the 
factorisation scale is found to move the peak of the Pt distribution about 0.2 GeV in the 
direction of the CDF data, with a slightly larger normalisation factor of 1.24. 

8 BR(W -> ev) = 0.1072 and BR(Z -> ee) = 0.03363 [22]. 
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Figure 4: distribution of Z bosons compared to CDF data [24 . 
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Figure 5: Pt distribution of Z bosons compared to D0 data |25j . 
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pp^ ev 




P T (GeV) 

Figure 6: Pt distribution of W bosons compared to D0 data [26J. 
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The small residual discrepancy between the solid lines in Figs. |H and |H1 and the data is 
easily understood. Note that the MRST2001 LO PDFs jTTj have been determined by a global fit 
to data using the conventional collinear approximation. A more precise treatment would fit the 
integrated PDFs, used as input to the last evolution step, to the proton structure function F 2 , 
for example, using the (z, /^-factorisation formalism at LO. We would expect this treatment 
to give slightly larger integrated PDFs than the conventional sets by a factor of ~ 1.1 and so 
eliminate the small discrepancy between the (z, /cj-factorisation predictions and the data. 

Alternatively, it was found in [§| that the major higher order corrections to the inclusive jet 
cross section in DIS could be accounted for by adding extra parton emissions to the LO diagram, 
7*g* — > q. In an axial gluon gauge, ladder-type diagrams gave the dominant contributions. 
Thus, it is likely that calculating the 0(as) subprocesses q* q 2 — > g V and q* g* — > qV using the 
(z, /^-factorisation prescription would account for any significant higher order corrections not 
already included and so reduce the observed discrepancy without refitting the input integrated 
PDFs. 

These reasons for the small discrepancies discussed in the previous two paragraphs can be 
illustrated by considering the proton structure function F 2 (x B , Q 2 ), which is plotted for two Q 2 
values in Fig. [7| In the collinear approximation, the LO prediction for this observable comes 
from 7*g — > q 

F 2 (x B ,Q 2 ) =^2e 2 q x B q(x B ,Q 2 ), (36) 
i 

indicated by the solid lines in Fig. [?[ which gives a good description of the data points since this 
data set was included in the MRST2001 LO global fit. The LO (z, kt) -factorisation prediction 
comes from 7*g* — > q and may be obtained from (85) of jU]: 

F 2 (x B ,Q 2 ) =^2e 2 g x B q(x B ,fi 2 )T q (fil,Q 2 ) + [ dz [ %B I X e 2 g f q (x, z, k 2 , fi 2 ), 

Jx Jul K t 1 X B p/X 

(37) 

where the Sudakov variables x = x + and j3 are given in (54) of and the factorisation scale \i 
is given in (25) of [Oj. The predictions of this formula are shown as the dashed lines in Fig. 
while the first term of (|37jl . representing the non-perturbative contributions from k t < fio, is 
also shown separately as the dotted lines. There is a clear difference between the predictions 
of (}3l)j) and ([37)1 . which increases as x B decreases, due to the extra kinematic factor in the 
second term of (|37|). This difference would be eliminated by fitting the input integrated PDFs 
using (JHZJ). Alternatively, a "NLO" prediction for F 2 may be calculated from the subprocesses 
1*9* ~^ QQ an d l*q* — > qg, and can be obtained from (94) of jHj- Here, a lower limit of is 
taken in the k' t integration, and the k' t < /i contribution is instead accounted for using the first 
term of (|37|). It is seen that these "NLO" predictions, shown as the dot-dashed lines in Fig. [7[ 
give almost the same results as (}36|) . 

5.2 Standard Model Higgs boson production at the LHC 

The Pt distribution (j3~Uj) for SM Higgs bosons of mass 125 GeV produced at the LHC (-y/i = 14 
TeV) is shown in Fig. [HI To allow direct comparison with the results of [30J, we do not account 
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Figure 7: Predictions for F 2 (x B ,Q 2 ) at Q 2 values of 22 GeV 2 (top) and 200 GeV 2 (bottom) 
in the collinear approximation, where the 1996/97 ZEUS data [2H| has been included in the 
MRST2001 LO fit 1 1J , and in the (z, fct)-factorisation approach using the same PDFs as input. 
The discrepancy between the data and the LO (z, k t ) -factorisation prediction can be eliminated 
by either refitting the input integrated PDFs or by adding some "NLO" contribution. 
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Figure 8: Pt distribution of SM Higgs bosons produced at the LHC with mass 125 GeV. 

for top quark mass effects. Note that the peak in the Higgs Pt distribution is broader and 
occurs at a higher Pt than for vector boson production. This is primarily due to the enhanced 
g — > gg colour factor (Ca = 3) compared to the q — ► qg colour factor (Cp = 4/3), resulting 
in a larger Sudakov suppression at low Pp. By the same reason the if-factor (J35|) is larger. 
For Pt < M^, the Pt distribution is in good agreement with recent, more sophisticated, 
resummation predictions (see, for example, [30J ) , bearing in mind the spread in the various 
predictions available due to the different approaches and PDFs used. However, the peak occurs 
at a Pt about 1-2 GeV lower than the majority of the resummation predictions (cf. Fig. 1 of 
[3*U]). Evaluating the total cross section by integrating over all Pt gives 38.6 pb, close to the 
next-to-next-to-leading order (NNLO) calculation which gives 39.4 pb [3*0] . 

Note that matrix-element corrections are necessary in parton shower simulations at large 
Pt- Without such corrections, the herwig parton shower prediction falls off dramatically at 
large Pt > Mh [Hill EH] ■ The fact that our prediction is much closer to the fixed-order results 
than parton shower predictions suggests that we have successfully accounted for a large part of 
the sub-leading terms. 
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6 Conclusions 



We have extended the method of (z, /^-factorisation using doubly-unintegrated parton dis- 
tributions [Oj to hadron-hadron collisions. The key idea is that the incoming partons to the 
subprocess have finite transverse momenta, which can be observed in the particles produced 
in the final state. This transverse momentum is generated perturbatively in the last evolution 
step, with a suitable extrapolation for the non-perturbative contribution. Virtual terms in the 
DGLAP equation are resummed into Sudakov form factors and angular-ordering constraints 
are applied which regulate soft gluon emission. By accounting for the precise kinematics in the 
subprocess, together with these Sudakov form factors and angular-ordering constraints, we are 
able to include the main part of conventional higher order calculations. 

We used this framework to calculate the Pt distributions of W and Z bosons produced at 
the Tevatron Run 1. The predictions gave a very good description of CDF and D0 data over the 
whole Pt range, after multiplying by an overall factor of 1.1-1.2, corresponding to multiplying 
each DUPDF by a factor < 1.1. We explained the origin of the need for this extra factor, which 
should not be regarded as a deficiency of our approach, but rather reflects the fact that the 
input integrated PDFs should themselves be determined from data using (z, /^-factorisation. 

We also used the framework to calculate the Pt distribution for SM Higgs bosons of mass 
125 GeV produced at the LHC. For P T < M H , our simple prescription was found to reproduce, 
to a fair degree, the predictions of more elaborate theoretical studies [3T)| . 
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